1 



o 
o 

(N 



X 

5-H 



On the width of collective excitations in chiral soliton models 

Herbert Weigel*) 
Fachbereich Physik, Siegen University, D-57068 Siegen, Germany 



In chiral soliton models for baryons the computation of hadronic decay widths of baryon 
resonances is a long standing problem. For the three flavor Skyrme model I present a solution 
to this problem that satisfies large-Afc consistency conditions. As an application I focus on 
the hadronic decay of the & and 0* pentaquarks. 

<^ \ §1. Statement of the problem 

OO . 

Hadronic decays of baryon resonances are commonly described by a Yukawa 
interaction of the generic structure 

Q-f Ant ~ 9i>B' <t>i>B , (l-i) 

p ; 

where B is the resonance that decays into baryon B and meson (j) and g is a coupling 
4h ■ constant. It is crucial that this interaction Lagrangian is linear in the meson field. 

If cf) is a pseudoscalar meson this interaction yields the decay width r{B' — > B(/>) oc 
ff 2 |p</>| 3 ] with p*0 being the momentum of the outgoing meson. 

The situation is significantly different in soliton models that are based on action 
functionals of only meson degrees of freedom, r = r[4>]. These action functionals 
contain classical (static) soliton solutions, ^ so i, that are identified as baryons. The 
interaction of these baryons with mesons is described by the (small) meson fluctua- 
tions about the soliton: <P = <P so \ + <fr. By pure definition we have 

^ ■ srm 



<5# 



= 0. (1-2) 



Thus there is no term linear in (p to be associated with the Yukawa interaction, 
eq. (1-1). This puzzle has become famous as the Yukawa problem in soliton mod- 
els. However, this does not mean that soliton models cannot describe resonance 
widths. On the contrary, these widths can be extracted from meson baryon scat- 
tering amplitudes, just as in experiment. In soliton models two-meson processes 
acquire contributions from the second order term 



r 



2 r 8 2 $ 



(1 ' 3) 



This expansion simultaneously represents an expansion in Nc, the number of color 
degrees of freedom: r = 0{N C ) and = O(N^). Terms 0(<f> 3 ) vanish in the limit 
Nq — ► oo. This implies that contains all large- Nc information about hadronic 
decays of resonances. We may reverse this statement to argue about any computation 
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of hadronic decay widths in soliton models: For Nc — ► oo it must identically match 
the result obtained from r^ 2 \ Unfortunately, the most prominent baryon resonance, 
the A isobar, becomes degenerate with the nucleon as Nc — ► oo. It is stable in that 
limit and its decay is not subject to the above described litmus-test. The situation 
is more interesting in soliton models for flavor SU(3). In the so-called rigid rotator 
approach (RRA), that generates baryon states as (flavor) rotational excitations of the 
soliton, exotic resonances emerge that dwell in the anti-decuplet representation of 
flavor SU(3). 1 ^ The most discussed (and disputed) such state is the + pentaquark 
with zero isospin and strangeness S = +1. In the limit Nc — > oo the (would-be) 
anti-decuplet states maintain a non-zero mass difference with respect to the nucleon. 
Therefore the properties of + predicted from any model treatment must also be seen 
in the 5-matrix for koan-nucleon scattering as computed from r^ 2 \ This (seemingly 
alternative) quantization of strangeness degrees of freedom is called the bound state 
approach (BSA) because in the S = — 1 sector the resulting equations of motion for 
<j) yield a bound state. Its occupation serves to describe the ordinary hyperons, A, 
U, £*, etc.. The above discussed litmus-test requires that the BSA and RRA give 
identical results for the <9 + properties as Nc — ► oo. This did not seem to be true 
and it was argued that the prediction of pentaquarks would be a mere artifact of the 
RRA. 2 ) Here we will show that this is a premature conclusion and that pentaquark 
states do indeed emerge in both approaches. Furthermore the comparison between 
the BSA and RRA provides an unambiguous computation of pentaquark widths: It 
differs substantially from previous approaches 3 ) that adopted transition operators 
for + — > KN from the axial current. 

This presentation is based on ref. 4 ) which the interested reader may want to 
consult for further details. 

§2. The model 

For simplicity we consider the Skyrme model 5 ) as a particular example for chiral 
soliton models. However, we stress that our qualitative results do indeed generalize 
to all chiral soliton models because these results solely reflect the treatment of the 
model degrees of freedom. 

Chiral soliton models are functionals of the chiral field, U, the non-linear real- 
ization of the pseudoscalar mesons*), <j) a 



U(x ,t) = exp 



— <f>a(x,t)\ a 
Jit 



(2-1) 



with A a being the Gell-Mann matrices of SU(S). For a convenient presentation of 
the model we split the action into three pieces 

r = r SK + r wz + r SB . (2-2) 



*-* A remark on notation: In what follows we adopt the convention that repeated indices are 
summed over in the range a, b, c, . . . = 1, . . . , 8, a, /3, 7, . . . = 4, . . . , 7 and k, . . . = 1, 2, 3. 
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The first term represents the Skyrme model action 



rsK 



+ 



1 



32e 2 



(2-3) 



Here /„■ = 93MeV is the pion decay constant and e is the dimensionless Skyrme 
parameter. In principle this is a free model parameter. The two-flavor version of 
the Skyrme model suggests to put e = 4.25 from reproducing the Z\-nucleon mass 
difference*^. The QCD anomaly is incorporated via the Wess-Zumino action 7 ) 

r W z = / <? x e^ P<TT tr [a^a u a p a a a T } , (2-4) 

with a M = U^d^U. Note that fwz vanishes in the two-flavor version of the model. 



The flavor symmetry breaking terms are contained in Isb 
f 2 

J 7T 



rsB = 



d xtr 



M ( U + C/t - 2 



M = diag (m 2 , m 2 , 2m 2 K 



m 



I) ■ (2-5) 



We do not include terms that distinguish between pion and kaon decay constants 
even though they differ by about 20% empirically. This omission is a matter of 
convenience and leads to an underestimation of symmetry breaking effects 8 ) which 
approximately can be accounted for by rescaling the kaon mass mx ^kSk/ fir- 

The action, eq. (2-2) allows for a topologically non-trivial classical solution, the 
famous hedgehog soliton 



^soi ~ U (x) = exp iX ■ xF(r) 



(2-6) 



embedded in the isospin subspace of flavor SU(3). The chiral angle, F(r) solves the 
classical equation of motion subject to the boundary condition F(0) — F(oo) = ir 
ensuring unit winding (baryon) number. The soliton can be constructed as a function 
of the dimensionless variable ef n r and is thus not subject to Nq scaling. 

In the RRA baryon states are generated by canonically quantizing collective 
coordinates A G SU (3) that describe the (spin) flavor orientation of the soliton, 
A(t)Uo(x)A* (t). The resultant eigenstates may be classified according to SU(3) 
multiplets; see ref. 9 ) for a review. 

§3. Large Nc P wave channel phase shifts with strangeness 

As motivated after eq. (1-3) we introduce fluctuations <f> ~ n a (x,t) 



U (x,t) = \/Uq{x) exp 



u 



X a r) a (x,t) 



(3-1) 



for the kaon fields. 10 ) Expanding the action in powers of these fluctuations is an 
expansion in rj a / f^ and thus a systematic series in 1/y/Nc- The term quadratic 



*^ To ensure that the (perturbative) n-point functions scale as n / 26 ) we substitute = 
93MeVVjVc73 and e = 4.25^/3/iVc in the study of the N c dependence. 
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in rj a describes meson scattering off a potential generated by the classical soliton, 
eq. (2-6). The P-wave mode is characterized by a single radial function 

In future we will omit the subscript that indicates the Fourier frequency. Upon 
quantization the components of the two-component iso-spinor x(u;) are elevated to 
creation- and annihilation operators. It is straightforward to deduce the Schrodinger 
type equation 

d? 2d 

h 2 r]{r) +uj[2\{r) - ujM K (r)]r](r) =0 with ft 2 = -_--_ + ^(r) . (3-3) 

The radial functions arise from the chiral angle F(r) and may be readily taken 
from the literature. 10 ) The equation of motion (3-3) is not invariant under particle 
conjugation u <-> — uj, and thus different for kaons (uj > 0) and anti-kaons (uj < 0). 
This difference stems from the Wess-Zumino term. Equation (3-3) has a bound state 
solution at u = —uja, that gives the mass difference between the /1-hyperon and the 
nucleon in the large-TVc limit. As this energy eigenvalue is negative it corresponds 
to a kaon, i.e. it carries strangeness S = —1. In the symmetric case (vtik = run) the 
bound state is simply the zero mode of SU(3) flavor symmetry. The WZ-term moves 
the potential bound state with S = +1 to the positive continuum and we expect a 
resonance structure in that channel. The corresponding phase shift is shown in the 
left panel of figure 1. No clear resonance structure is visible; the phase shifts hardly 
reach it/ 2. The absence of such a resonance has previously lead to the premature 
criticism that there would not exist a bound pentaquark in the large-A^c limit. 2 ) 



§4. Constraint fluctuations 



To study the coupling between the fluctuations and the collective excitations we 
generalize eq. (3-1) to 



U(x,t) = A(t) v / ?7o(^)exp 



—X a r] a (x,t) 

Jn 



(4-1) 



BOUND STATE APPROACH 
m K = 495 MeV 



Fig. 1. Large TVc P-wave phase shifts with strangeness S = +1 as function of the kaon momentum. 
Left panel: unconstrained; right panel: constrained to be orthogonal to the collective rotation. 
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These fluctuations dwell in the intrinsic system as they rotate along with the soliton. 
The kaon P-wave is subject to the modified integro-differential equation 

" POO 

h 2 rj(r) + oo [2A(r) - ooM K (r)\ rj(r) = -z(r) / r' 2 dr' z(r')2X(r')7j{r') 
x 2A(r) - (oo + oo )M K (r) - oo 



+ 3.) (2A(r) - oo M K (r)) 

000 — 00 00 J 



, (4-2) 



for the flavor symmetric case*). The radial function n(r) is defined according to 
eq. (3-2) and z(r) = 



■ sin 



F(r) 



is the collective mode wave-function nor- 
malized with respect to the moment of inertia for flavor rotations into strangeness 
direction, Ok = f 2 J d 3 rMx{r) sin 2 = 0(Nc)- The non-local terms without 
Xj\,e reflect the constraint J drr 2 z(r)Mx (r)rj(r) = which avoids double counting 
of rotational modes in strangeness direction. The interesting coupling is contained 
in the interaction Hamiltonian 



H 



int 



d 



K 



ia/3 D^aRp 



I 



d 3 r z(r) [2A(r) - oo M K (r)} x ^(x, t) , (4-3) 



where £ a = D a 6% are the fluctuations in the laboratory frame, that we actually 
detect in KN scattering. The collective coordinates are parameterized via the adjoint 
representation D a b(A) = \ tr [Aq^A;,^] and the SU(3) generators R a . Integrating 
out the collective degrees of freedom by means of standard perturbation theory 
induces the separable potential 



K0|F int |(iTAO /=o >| 2 , \{A\H int \(KN) I=0 



000 



+ 



00 



00 A + CO 



(4-4) 



These matrix elements concern the T-matrix elements in the laboratory frame. Since 
the laboratory and intrinsic T-matrix elements are identical for the + channel, 11 ) 
we may add the exchange potential, eq. (4-4) in the intrinsic frame. We define matrix 
elements of collective coordinate operators 



{e + \d Za pD +a R p \n) =:X 



N C , 

and 

32 



(A\d^D^ a R p \p) 



(4-5) 



to end up with eq. (4-2). The first factor in the coefficient ooq = 2 



N C 
A0 h 




arises in the equation of motion because the potential, eq. (4-4) is quadratic in 
the fluctuations. The remaining (squared) factors stem from the definitions of Xq^a 
and the constant of proportionality in H[ nt . The Xg^ must be computed with the 
methods provided in ref. 12 ) but generalized to arbitrary (odd) Nc-^ For Nc — ► oo 
we have X0 — > 1 and Xa — > 0. From the orthogonality conditions of the equation of 
motion (3-3) we straightforwardly verify that 



rj(r) = rj(r) — az(r) 



with 



drr 2 z(r)MK{r)rj(r) 



(4-6) 



The more complicated case uik 7^ is at length discussed in ref. 4 - 1 
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solves eq. (4-2) for large Nc- This is essential because, as z(r) is localized in space, 
r\ and rj have identical phase shifts! Hence the litmus-test discussed in the introduc- 
tion is indeed satisfied. The physics behind rj is best understood when introducing 
reduced wave-functions rj(r) that solve eq. (4-2) modified with Xq = Xy\ = 0. i. e. 
the collective excitations are decoupled. These wave-functions are still orthogonal 
to the collective modes and actually lead to the background phase shift shown in 
figure 1. Having obtained the rj(r) we may again switch on the exchange contribu- 
tions, eq. (4-4). The additional separable potential is treated by standard R-matrix 
techniques and augments the phase shift by 



tan(5 R (fc)) 



r(w fe )/2 



lvq - Uk + A(uj k ) 



(4-7) 



^o^ 3 is the RRA result for the excitation energy of 0. This phase shift 



Here loq 

exhibits the canonical resonance structure with the width and the pole shift 



r{u k ) = 2ku X% 



I 

Jo 



r 2 dr z{r)2\(r)rj (r) 



o 

oo 



i r 00 

A(u k ) = - V / qdq 

2nuJ k J 



— LO q LO k + UJ q 



(4-8) 
(4-9) 



respectively. We have numerically verified that in the large-iVc limit with Xq = 1, 
the phase shift from eq. (4-7) is identical to what is labeled resonance phase shift in 
figure 1, that we calculated as the difference between the total {rj) and background (rj) 
phase shifts. For finite Nq we have Xq ^ 1 and X\ ^ so the R-matrix formalism 
becomes two-dimensional (A and + exchange). Contrary to earlier criticisms 2 ) the 
large Nc pentaquark phase shift indeed resonates! 

§5. Results 



In figure 2 we show the resonance phase shift computed from eq. (4-7) for various 
values of Nc- First we observe that the resonance position quickly moves towards 
larger energies as Nc decreases. This is mainly due to the strong Nc dependence 
of ujq: For Nc = 3 it is twice as large as in the limit Nc — ► oo. The pole shift A 
is actually quite small (some ten MeV) so that ujq is indeed a reliable estimate of 
the resonance energy. Second, the resonance becomes shaper as Nc decreases. To 
major parts this is caused by the reduction of Xq. 

We now turn to more quantitative results for which we also include flavor sym- 
metry breaking effects. Then the resonance position changes to 



LOQ 



_3T_ 
20^ 



+ w 



+ o 



1 



(5-1) 



where r = 0(Nq) is a functional of the soliton that is proportional to the meson 
mass difference, m^—m 2 . The O (l/Nc) piece is sizable for Nc = 3 and we compute 
it in the scenario of ref. 12 ) We then find ujq « 700MeV; taking model dependencies 
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into account we expect the pentaquark to be about 600 . . . 900MeV heavier than the 
nucleon. 

For the width calculations there are two principle differences. First, the inter- 
action Hamiltonian acquires an additional term 

tffnt = ( m K ~ m l) d ia pD ja D 8 p J d 3 r z{r)i{r%(x ,t) Xi , (5-2) 

The radial function j(r) is again given in terms of the chiral angle. 4 ) Second, the Xa 
does not vanish as Nq — ► oo and the i?-matrix formalism is always two dimensional. 
Nevertheless, the large-TV^* solution is always of the form (4-6) and the BSA phase 
shift is recovered. The width function turns to 

2 

(5-3) 

where Xq and Yq = y/8Nc /3(<9 + \d3 a pD +a D%p\n) are to be computed in the RRA 
approach with full inclusion of flavor symmetry breaking effects. 

This width function is shown (for Nq = 3) in figure 3 for and its isovector 
partner 0* . The latter merely requires the appropriate modification of the matrix 
elements in eq. (4-5). Most importantly, the k 3 behavior of the width function, as 
suggested by the model, eq. (IT) is reproduced only right above threshold, afterwards 
it levels off. Second, and somewhat surprising, the width of the non-ground state 
pentaquark is smaller than that of the lowest lying pentaquark. Our particular 
model yields 40MeV and i~e* k, 20MeV. We note that there are certainly 

model ambiguities in these results. 



r(uj k ) = 2ku; 



2 dr z(r) 



2X X(r) + ^ ( 



ml) 




Fig. 3. Model prediction for the width , r(u>) of (9+ (left) and 6>* + (right) for N c = 3 as function 
of the momentum k = \/co 2 — m 2 K for three values of the kaon mass. Note the unequal scales. 
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§6. Conclusions 



We have discussed the chiral soliton model approach to KN scattering in the 
S = +1 channel which contains the potential 14 ) + pentaquark, a state predicted as a 
flavor rotational excitation of the soliton. Though the exactly known large Nc phase 
shift suggests otherwise, the © emerges as a genuine resonance. Our central result 
is the width function for © — ► KN. In the flavor symmetric case it contains only a 
single collective coordinate operator and is thus very different from estimates that 
extract an effective Yukawa coupling from the axial current matrix element. 3 ) Since 
our approach matches the exact large Nc result, we must conclude that those axial 
current scenarios are erroneous 13 ) and that the cancellation among contributions to 
this matrix element is an invalid argument for a small pentaquark width. 
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